Accepted by ApJ Letters 

Preprint typeset using LAT^X style emulateapj v. 21/08/00 



A RELATIVISTIC JET IN THE RADIO-QUIET QUASAR PG 1407+263 
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ABSTRACT 

We present the results of a multi-epoch radio monitoring campaign measuring the milliarcsecond 
structure of the jet in the radio-quiet quasar PG 1407+263. This is the highest-sensitivity, highest- 
resolution multi-year study of a distant active galaxy. The observations are naturally explained in terms 
of a beamed relativistic jet, some of whose fluctuations in flux density can be ascribed to interaction 
with the narrow-line region of the quasar. The optical properties of PG 1407+263, in particular the low 
equivalent widths of the emission lines, may be related to the fact that we are viewing this quasar almost 
pole-on, giving us a direct view into its broad-line region. 

Subject headings: galaxies: jets — quasars: individual (PG1407+263) — radio continuum: galaxies 



1. INTRODUCTION 

In recent years, the weak radio core emission associated 
with the class of so-called 'radio-quiet' quasars (for exam- 
ple, Kellermann et al 1989; Miller et al 1990, 1993) has 
been shown by a succession of Very Long Baseline Inter- 
fcrometry (VLBI) imaging experiments usually to be com- 
pact on milliarcsecond (mas) scales (Blundell et al 1996; 
Falcke, Patnaik & Sherwood 1996; Blundell & Beasley 
1998a). The typical physical size of < a few cubic par- 
sees for these radio luminosities and corresponding bright- 
ness temperatures (> 10 6 K) strongly argues against star- 
formation as the origin of this radio emission, despite the 
striking continuity of the correlation of far-infra-red and 
radio emission from star-forming galaxies to radio-quiet 
quasars (Sopp and Alexander 1991a, b). Instead, a mech- 
anism resembling that which powers the core emission in 
the more radio-luminous 'radio-loud' quasars is suggested. 

One quasar in particular, PG 1407+263, identified in 
a mas imaging survey of radio-quiet quasars (Blundell 
& Beasley 1998a) clearly has multiple components and 
seemed an important target to monitor subsequently with 
VLBI techniques. It was first discovered as a quasar in 
the Palomar-Green survey (Schmidt & Green 1983) with 
z = 0.94 ± 0.02 (McDowell et al. 1995). With a Hubble 
constant of 70 kms -1 Mpc _1 , and (0 M ,fi A ) = (0.3,0.7), 
the luminosity distance ss 6.2 Gpc and lmas « 7.9 pc. 



2. MEASUREMENTS 

We observed PG 1407+263 at 8.4 GHz using the Very 
Long Baseline Array (VLBA) with the Very Large Array 
(VLA) in phased array mode on 1996 June 15 (VLBA 
only), 1997 April 27, 1998 January 31, 1999 March 20 
and 1999 July 20. We usually recorded a 64-MHz band- 
width (a 32-MHz bandwidth was used for epochs 4 & 5) 
with 2-bit sampling and dual circular polarization. The 
observations were phase-referenced (Beasley and Conway 
1995), i.e. regular observations of a bright extragalactic 
source close on the sky to the target were made to derive 
residual antenna-based and atmospheric phase corrections, 
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which were interpolated to calibrate the target data and 
provide accurate astrometric referencing. The calibrator 
source used was OQ 208, separated from PG 1407+263 by 
2.2°. We typically used a switching cycle of two minutes 
on the calibrator and three minutes on the target. The 
total on-source time for OQ208 was 0.5 hrs (epoch 1) and 
2.8 hrs (epochs 2, 3, 4 and 5). Throughout the three-year 
monitoring, a pair of < mJy components were persistently 
detected (Figure 1), albeit with changes in structure and 
flux density (Figure 2). 

3. LIGHT CURVES AND STRUCTURAL CHANGES 

One striking change seen between successive epochs is 
that the structure of each component appears to vary: this 
is particularly marked in the final three epochs (Figure 1). 
Real changes in shape represent an intriguing aspect of 
the parsec-scale evolution of the source: since 1 mas cor- 
responds to 7.9 pc, real expansion between epochs of one 
tenth of a synthesized beam would correspond to appar- 
ent motion of a few times the speed of light. The com- 
ponents appear both resolved and unresolved at different 
epochs. Component B appears to be significantly resolved 
(with a high total flux density) in earlier epochs, while 
both components appear essentially unresolved (with low 
flux densities) in the final epoch. Unmodeled residual spa- 
tial gradients in tropospheric (and ionospheric) delays lead 
to interferometer phase errors which are cancelled towards 
the bright calibrator source OQ 208, but remain in the tar- 
get data. Thus, apart from the above general assertions 
above about component B, it is difficult to draw solid con- 
clusions from the structural variations: simple centroid 
fitting to the components, while technically possible, can 
be ambiguous in interpretation. 

The rms variation of the position of component A is 
» 300/i-arcsec. Typical estimates of the effects of tropo- 
spheric delay errors on phase-referenced VLBI data with 
a calibrator 2.2 deg away are 100-200 /^-arcsec. In addi- 
tion, the typical signal-to-noise (S/N) of the detections 
is 5-10 to 1; this would introduce a positional variation 
roughly equal to the synthesized beam divided by the S /N 
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- i.e. 100 to 250 /z-arcsec. Therefore, the position wan- 
dering largely results from a combination of unmodeled 
tropospheric delay errors and S/N considerations. Given 
the weak (< mJy) flux density levels in these images, it 
is not possible to correct for these residual errors. During 
the course of these observations the relative position shift 
of components A and B appeared at times to be system- 
atic, possibly indicating superluminal motion (Blundell & 
Beasley 1998b). Analysis of the existing data set does not 
allow a definitive determination of superluminal motion in 
this radio-quiet quasar, although we in any case infer it to 
have a highly relativistic jet speed (Doppler factor >10) in 
the analysis below. The rapid variations in component B's 
flux density are unlikely to result from imaging coherence 
loss, which would affect both components; A's relatively 
constant peak flux density implies that B's flux density 
variations are intrinsic. The dramatic changes in B's flux 
density occur on timescales of 100-300 days, reminiscent 
of those seen in the cores of radio-loud quasars (Hough et 
al 1999). 
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Fig. 1.— VLBA images of PG 1407+263 restored with the same 
synthesized beam, 1.8 X 0.9 mas with position angle 170 deg East 
of North. The contours in each case arc ±0.3, 0.42, 0.6, 0.84, 1.2, 
1.68, 2.4mJy/beam. 



4. RECENT RADIO-LOUDNESS OF PG 1407+263 

In order to establish whether PG 1407+263 is truly a 
radio-quiet quasar, we made deep VLA images to search 
for evidence of extended radio structure or hotspots. 
These images revealed only independent sources including 
an apparent double 53 arcsec east of the quasar. Keller- 
mann ct al (1989) commented that this double structure 
might be related to the output of the quasar, since it lacks 
any optical identification on the Palomar Observatory Sky 
Survey. We made a more sensitive test by observing the 
field in the near infra-red fC-band on the UK Infra-Red 
Telescope (UKIRT). Figure 3 indicates each part of the 
'double' radio structure is identified with a galaxy and is 
thus unrelated to the quasar. Therefore, we do not be- 
lieve that PG 1407+263 has been radio-loud for the last 
10 6-7 years, where this timescale is taken from the typical 
radiative lifetimes for synchrotron-emitting particles in the 
lobes of radio-loud quasars (Blundell & Rawlings 2000). 
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Fig. 2. — Peak flux densities at 8.4 GHz of components A & B 
at each epoch plotted against time. The dotted line corresponds 
to the eastern component (A) and the solid line corresponds to the 
western component (B). 

5. PHYSICAL INTERPRETATION 

It is feasible that the observed structural changes in 
component B are related to at least one unresolved sub- 
component. The much closer source 3C 120 (Gomez et 
al 2000) is a good example of how jet knot components 
can appear and fade. Gomez ct al have interpreted this 
in terms of interaction between the jet and the inhomo- 
geneous ambient medium. Observed at lower resolution, 
the jet in 3C 120 would present apparently structurally 
changing partially resolved components. If we are indeed 
observing such a phenomenon, then the offset between the 
two local surface brightness maxima in component B may 
give us some idea of the scale of the jet diameter. The per- 
pendicular offset between the centroid of the two B sub- 
components and the line connecting the two main compo- 
nents, A and B is about 0.7mas w 5.5pc. Therefore, in the 
following, we adopt a fiducial scale for the jet diameter, D 
of 5pc. The causal limit on the diameter deduced from the 
variability time scale, At, and corresponding to a Doppler 

factor 5 is: D < c(l + z)- 1 At 5* 0.4 ( 1?8 ^ s ) (&) pc. 

The variation from epochs 1 to 2 entails a factor ~ 2 de- 
crease in flux density over a timescale w 400days. This im- 
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plies an upper limit on the diameter w 1.6pc for a Doppler 
factor of 10. This is intriguingly less than, but of the same 
order as, the approximate scale size estimated from the 
VLBA images. Of course, the discrepancy decreases with 
increasing Doppler factor. [If the jet is pole-on so that 
S k 2T (where T is the bulk Lorentz factor), given the ex- 
istence of apparent proper motions ~ 10 c in some quasars, 
then 5 ~ 20 is not out of the question.] Clearly, the rec- 
onciliation of these two scales is in the direction of larger 
Doppler factors, not smaller. This is an indication that 
this source is highly beamed. 

A second, more persuasive, indication for beaming 
comes from estimating the intrinsic source power. The 
emitted power per steradian is: J j' v dV = (5~( 3+a ) (1 + 
^-(l-a) £)2 Hz -1 sr _1 , where the integral is that 
of the rest-frame emissivity, j' v , at the observed frequency 
v over the proper volume (V) of each blob and F v is the 
flux density at v. 




14 h 09 m 28.5 s 28.0 s 27.5 s 27.0 s 

Right Ascension (J2000) 

Fig. 3. — The top image is a seven arcmin-square VLA image at 
8.4 GHz in D-array. The lower image shows a grey-scale of near-IR 
K-h&rid on which are overlaid contours at 8.4 GHz from the VLA A- 
array, the lowest of which is 0.04 mjy/beam. The two radio sources 
are clearly identified with resolved 18.5-magnitude galaxies, indicat- 
ing they are independent radio structures from the quasar. 

For unity Doppler factor, the total emitted power, 
at 8.4 GHz, of both components at epoch 1 would be 
1.1 x 10 24 W Hz -1 sr _1 . This is comparable to the core 
powers of the most powerful galaxies in the B2 sample 
(de Ruiter et al 1990). If PG 1407+263 were similar to 
the B2 sources then, using the statistical relationship be- 
tween core and extended powers derived by de Ruiter et al 
(1990), its corresponding 1.4 GHz extended power would 
be of order 10 28 W Hz" 1 and the 8.4 GHz flux density 



would be of order one Jansky. However, since the flux 
density of any extended emission from PG 1407+263 is 
less than O.lmJy this is another indication that the core 
is beamed. The Doppler factor has to be increased to ~ 10 
in order to reduce the expected extended flux density to 
the observed upper limit. The additional parameters de- 
scribing each knot are the spectral index a, the electron 
energy index, a — 2a + 1, the volume, V, of the knot 
in its rest frame, the Doppler factor, 6 and the ratio, c\e, 
of energy in relativistic electrons and positrons to other 
particles. The rest-frame magnetic field that minimizes 
the total energy density of a spherical blob in the rest 

/ 1 — rv \ 

frame is: B'. = e 



*D 2 L V 



[(a + 1) C^ L (a) (1 + c E ) m~ L c /(a, 7l , 72 ) F v v a 

Here Ci{a) is a coefficient involving T-functions that ap- 
pears in the expression for the angle-averaged synchrotron 
emissivity; e, m e and c are the electronic charge, electronic 
mass and speed of light respectively, and /(a, 71, 72) = 
(a - 2)- 1 7 2 ~ Q (l - (7 2 /7i) 2 ~ a ) 4 - We have estimated 
the minimum-energy magnetic field and corresponding 
particle energy densities, from the flux density of com- 
ponent B at epoch 1 and for proper jet diameters of 
0.375, 0.75 and 1.5 mas as functions of Doppler factor. 
Figure 4 shows the associated jet energy flux, Fe ~ 

where A\ et is the 



C/Mjct, w licit; /ij e t 

cross-sectional area of the jet, p is the rest-mass den- 
sity, and where the Poynting flux of a perpendicular mag- 
netic field has been included. We assume that the jet is 
pole-on and purely relativistic, so that pc 2 /Ap <C 1 and 
ce — 0. The almost constant value of the jet energy flux 
for Doppler factors > 2 is the result of the (5 _2 -dependence 
of the magnetic and particle energy densities canceling 
the r 2 -dependence of the energy flux. If we assume that 
0.75 mas is a reasonable upper limit to the jet diameter, 
then the jet is no more powerful than about 7 x 10 36 W. 
For jet diameters less than the fiducial value the jet is less 
powerful, albeit still relativistic. We note that this is not 
the first radio-weak object for which relativistic velocities 
have been inferred. Brunthalcr et al. (2000) inferred su- 
perluminal motion from their observations of the Seyfert 
galaxy IIIZw2, directly implying relativistic motion. 
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Fig. 4. — Jet energy flux as a function of Doppler factor for three 
assumed values of the jet diameter: the fiducial jet diameter and a 
factor of two above and below this. 
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6. DISCUSSION 

The low extended emission relative to the core in 
PG 1407+263 is a good argument that this jet is relativis- 
tically beamed. If our inference of partially resolved struc- 
ture in the jet is correct then the case for relativistic beam- 
ing is reinforced by the variability time scale (~ 400 days 
for epochs 1-2). If the jet diameter is indeed ~ 5 pc then, 
given the typical expansion rates of jets, the components 
are of the order of 100 pc from the core, placing them in 
the narrow-line region. These two arguments point to a 
Doppler factor of > 10. As Figure 4 shows, the dependence 
of energy flux on Doppler factor and jet radius indicates 
a kinetic power no greater than about 10 37 W. Using the- 
oretical (Bicknell et al. 1998) and observational (Willott 
et al 1999) relations between jet and radio power, such an 
energy flux would be typical of a radio galaxy jet no more 
than an order of magnitude above the FRI/FRII break 
(Fanaroff & Riley 1974). The interest that attends this 
galaxy is the existence of an FRI radio power associated 
with a bright quasar so it is unsurprising that this quasar 
is classified as radio quiet (Blundell & Rawlings 2001). It 
appears the only reason that we detect any radio flux is 
that it is one of the few PG quasars that is beamed towards 
us. 

As well as its variation between epochs 1 and 2, com- 
ponent B also exhibits a significant rise in flux density F v 
between epochs 3 and 4 then a sharp drop (At — 122 days) 
between epochs 4 and 5 with AF„/F U ~ 1. This variabil- 
ity indicates a diameter < 0.5(5/10) pc for component B 
around this epoch - much less than the indicative jet di- 
ameter ~ 5 pc inferred above. There are two possibilities: 
(1) the jet may not be as wide as 5 pc or (2) the variation is 



due to a small section of the jet being affected by its inter- 
action with clouds in the narrow-line region of the quasar. 
This second case is similar to what has been observed in 
3C 120 by Gomez et al (2000). Our observations therefore 
suggest a relativistic jet, of only moderate power, propa- 
gating through the clumpy narrow-line region of a bright 
host quasar. 

Docs our consideration of the radio structure of 
PG 1407+263 shed any light on its curious optical prop- 
erties e.g., the low equivalent width of the emission lines 
(McDowell et al. 1995)? If the Doppler factor of the jet is 
as high as 10, then we are viewing PG 1407+263 within a 
few degrees of pole-on and looking directly into its broad- 
line region. Thus, the emission line spectrum may be di- 
luted by a direct line of sight to the accretion disk con- 
tinuum. An argument against this is that the PG quasars 
show no evidence for orientation-dependent effects and this 
is generally ascribed to the initial selection via ultra-violet 
excess. However, supporting evidence for PG 1407+263 
being close to pole-on is the large velocity dispersion of its 
emission line gas w 10,000kms _1 (McDowell et al. 1995), 
indicating that we may be looking directly into the inner 
part of the broad-line region. 
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